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Coodness-of-£it  tests  are  developed  for  the  gamma  dis¬ 
tribution  when  the  scale  and  location  parameters  are  unspec- 
i l i o j  ant  must  he  estimated  from  the  sample  data. ^The 
Kol  mo  jorov-Sni  rnov,  Cramer-von  /lises,  and  Anderson-Darling 
statistics  are  used  to  develop  the  tables.  a  comprehensive 
power  study  is  conducted  to  compare  the  Kolmogorov-Smirnov, 
Cramer-von  mises,  and  Anderson-Darling  goodness-o f - f i t 
tests.  An  analysis  is  performed  to  determine  the  functional 
rel<*  t  ionshi  ,>  between  the  critical  values  and  the  shape 
parameter  . 
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Abstract 


Vhe  ftnderson-Darling,  Craner-von  ,'dses,  and  the 
Kol : logorov-Sn  i  r  nov  statistics  are  used  to  develop  a  new  test 
oi  fit  for  trie  tiiree-.:ara:;eter  'jam me  distribution  with 
unknown  silage  and  location  parameters.  The  critical  values 
generated  were  obtained  by  a  Monte  Carlo  procedure,  tor 
earn  value  ot  n  ( s  a  m  p  1  e  size),  saw  pie*  sets  were  drawn 

fro w  a  g  < . i a  population  wnosa  shape  is  specified.  The 

location  ana  scale  parameters  are  estimated  from  the  data, 
and  t;.;-.  turec  statistics  are  calculated  based  on  the  esti¬ 
mated  distribution.  The  simulation  was  preformed  for  sample 

r>  *3  •  >.3 

sizes  n  =  5#  Tt* ,  15.. .3(4- and  shape  parameters,  K  =  .5,  byt'-, 

1 . 5  . . .  4  . 

Using  gamma  distributions  for  shape  coual  to  1.5  and 

U, 

•Svft,  the  power  of  each  test  is  investigated  against  ten 

alternative  distributions  for  sample  sizes  n  =  5,  15,  and 

3  ■ 

30.  In  general  both  the  Anderson-Darling  and  the  Craner-von 
ibises  tests  are  more  powerful  than  the  Kolnogorov-Smirnov 
test.  Except  for  the  case  where  the  alternative  distribu¬ 
tion  is  lognormal,  the  Craner-von  Pises  test  is  the  most 
powerful  test. 

The  functional  relationship  between  tno  critical  values 
of  the  Anderson-Darling,  Crow.er-von  mises,  and  Kol mogorov- 

Siairnov  is  also  examined.  A  critical  value  for  a  shape 

U .  c 

parameter  between  1.5  and  -4rrt  which  is  not  included  in  the 
tables  can  tnen  be  easily  derived  from  this  functional 

relationship. 
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A  MODIFIED  KOLiiGCOKOV-SMIKNOV,  CRAMER— 
VO N  WISES  AND  ANDERSON-DARLING  TEST 
FOR  THE  GAMMA  DISTKIUUTION  tvITl. 
UNKNOWN  LOCATION  AND  SCALE 


I .  In  trod uc  t  ion 

Currently  ttie  U.S.  Air  Force  is  placing  more  and  more 
emphasis  on  system  availability,  maintainability,  and  reli¬ 
ability,  both  in  research  and  development  and  in  day  to  day 
operations.  Of  particular  importance  to  the  Air  Force  is 
the  ability  to  predict  ti me-to-fai lure  of  equipment. 
Studies  in  probability  and  statistics  have  increased  under¬ 
standing  of  some  Key  probability  distributions  used  in  pre¬ 
dictin';;  t i me- to- fa i lure*  Among  the  most  commonly  used 
continuous  distr  ibutions  in  this  area  are  the  beta,  gamma, 
exponential,  Wei  bull,  and  lognormal  distributions. 

Often  in  these  studies,  analysts  are  confronted  with 
the  problem  of  testing  agreement  between  probability  theory 
and  actual  observations.  In  other  words,  given  n  observa¬ 
tions  of  some  variable,  say  time-to-f ailure,  the  problem  is 
to  find  out  if  it  can  be  regarded  as  a  random  variable 
bavin j  a  given  probability  distribution.  The  general 
approach  to  tiie  solution  of  this  problem  is  known  as  the 
goodness-oL-t  it  test.  In  more  precise  terms  let  xj/*2**,xn 
be  a  random  sample.  Then  a  statement  of  the  joodness-of-£ i t 


I‘0:  1  (x)  =  lH(x) 


Ha:  F(x)  /  F„(x) 

where  i-'(x)  is  the  actual  distribution  function  of 
iH(x)  is  the  hypothesized  distribution  function. 


VF 


and 
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Two  coinuonly  used  ijoodness-of-f  i  t  teats  are  the  Chi- 
square  test  and  the  Kol moyo rov-Smi rnov  test.  The  Chi-square 
test  compares  observed  frequencies  with  expected  frequencies 
oi  tin!  hypothesized  distribution.  It  is  restricted  to  larqe 
sail [  lcs — cpi  rox i matel y  25  or  qreater  (2:73).  The  Kolnoqorov- 
hi.tirnov  (h-s)  test  compares  cumulative  frequencies  between 
the  actual  sample  usinq  a  step  function,  aqainst  correspond- 
ia  :  values  usinq  tut  hypothesized  cumulative  distribution 
function.  The  K-S  test  can  be  used  for  larqe  or  small 
samples;  however,  it  is  restricted  to  distributions  which 
are  lolly  specified.  li.  w.  Lilliefors  developed  a  rjoodness- 
ot-fit  test  for  the  normal  (19),  and  exponential  distribu¬ 
tion  (22/),  which  can  be  used  for  small  samples  where  the 
k.a  r  a  i>.  e  t  e  r  s  must  be  estimated  from  the  sample  data.  When 
parameters  are  estimated  from  sample  data,  the  test  is  said 
to  be  a  modified  test. 

followinq  Lilliefors'  technique,  several  other  modified 
tests  have  been  documented.  U.  Cortes  developed  a  modified 
Kolo^orov-bmirnov  test  for  the  three-parameter  Weibull  and 
qamma  distributions  (5).  J.  Bush  expanded  the  qoodness-of- 
lit  test  for  the  Weibull  to  include  the  modified  Craner-von 
wises  (.  '•)  and  anderson-bar  1  inq  (a2)  te^ts  (3h).  The  modi- 
Lien  . of moqorov-  'Jmirnov,  Cramer-von  lises,  unu  Anderson- 
barlinq  tests  liave  also  been  clone  for  the  uniform,  normal, 
uapi.iCe,  ex.onentiul  and  Caucuy  distributions  (11).  In  1973 
.•I  a  an,  bc.ieuer,  and  fertiq  dcsijned  two  new  test  statistics 
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called  the  L  and  S  statistics.  The  L  and  S  test  statistics 
were  used  to  develop  a  goodness-of-fit  test  for  the  two 
parameter  Wcibull  with  unknown  parameters  (22). 

In  lr;Cl  Koutrouvelis  and  Kellerr.ieier  introduced  a 
goodncss-of-fit  test  based  on  the  empirical  chacteristic 
function  when  the  parameters  must  be  estimated  (17).  This 
test  statistic  could  be  used  as  an  alternative  to  the  f;DF 
statistic  if  the  characteristic  function  is  more'  easily 
determined  than  the  distribution  function. 
limp i  r  ical  ui  str  ibution  Function  Statistics 

A  general  class  of  statistics  used  for  the  goodness-of- 
fit  tests  is  called  empirical  distribution  function  (EDF) 
statistic^.  historically  SDI  statistics  have  been  used  in 
cases  where  the  parameters  are  either  known  or  unknown.  In 
most  instances  EDF  statistics  are  easily  calculated  and  are 
competitive  in  terms  of  power.  This  class  of  statistics  is 
based  on  a  comparison  between  the  cumulative  distribution 
function,  J  (x),  and  the  empirical  cumulative  distribution 


function  S  (x)  defined  as 


no.  of  x^  <  x 


The  test  procedure  is  summarized  as  follows:  given  a  sample 
from  some  population,  the  £DF  tests  reject  !i(J:F(x)  =  F,,(x) 
when  the  difference  between  F  (x)  and  Sn(x)  is  large.  here 
Fjt(x)  is  the  hypothesized  distribution.  In  general,  LDF 
tests  are  valid  when  the  distribution  is  fully  specified, 
iiowever,  bavid  and  Johnson  showed  that  the  oistrioution  of 
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an  EbF  statistic  depends  only  on  the  functional  form  of  the 
distribution  and  not  on  the  unknown  parameters  when  the 
estimated  parameters  are  location  and  scale  (6).  It  is  this 
principle  that  permits  us  to  generate  valid  critical  value 
tables  for  the  gamma  distribution  which  depend  only  on  the 
shape  parameter  and  sample  size. 

It  is  important  to  note  that  this  thesis  uses  a  modi¬ 
fied  form  of  the  EOF  statistic,  because  the  cumulative 
distribution  is  not  fully  specified.  An  estimated  distribu¬ 
tion  function  is  used  whose  parameters  are  derived  from  the 
observed  sample. 

The  Kolmogorov-Sni  rnov  Statistic 

The  KolDoqorov-Srairnov  (K-S)  statistic  is  defined  as 
the  absolute  value  of  the  difference  between  F(x)  and  Sn(x) 
or , 

D  =  I  F (x)  -  Sn(x) I .  (2) 
In  using  the  K-S  statistic  for  the  goodness-of-f it  test,  we 
are  interested  in  the  greatest  absolute  difference  between 
I  (x)  and  bfi  ( x )  (?) .  Therefore  the  test  statistic  is 

T  *  syp| F(x)-Sn (x) | .  (3) 
The  Anderson -Dari ing  Statistic 

It  is  known  that  goodness-of-f it  tests  which  use  actual 
observations  without  grouping  are  sensitive  to  discrepancies 
at  the  tails  of  the  distribution  rather  than  near  the  median 
(?f>:2).  The  Anderson-Darling  test  statistic  overcomes  this 
problem  by  accentuating  the  values  of  Sn(x)  -  F(x)  where  the 
test  statistic  is  desired  to  have  sensitivity.  More 


IT* 


specifically,  tne  Anderson-Darling  statistic  is  based  on  a 
weighted  average  of  the  squared  discrepancy,  (i.e.  [ S n  ( x )  - 
F(x)]2  weighted  by  ♦( F(x))  or 

An2  *  ny**[Sn(x)-F(x)  J2  *  (F(x))dP(x)  ,  (4) 

where 

♦  ( i’  ( x )  )  =  [  P  ( X )  *  (1-F  (x)  )  J"1.  (5) 


Using  the  computational  form 
n 

An2  -  -n  -  A  £(2j-l)[ln  F(xj)  +  In  (1-Fn„j+1)], 

n  j  =  1  J 


(6) 


the  test  procedure  is  as  follows: 

1)  Let  x1<x2<...<xn  be  n  observations  in  the  sample. 

2)  Compute  An2 

3)  If  An2  is  too  large,  the  hypothesis  is  to  be  rejected. 
The  Cr amc-r-von  Mises  Statistic 


The  Cramer-von  Mises  statistic  is  a  special  case  of  the 
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A(1  c  with 


[  F  ( x )  ]  =  1  and  is  written  as 

Wn2  =  n(x)  "  F(XH2  dx»  (7) 

"—CD 

This  test  procedure  is  the  same  as  outlined  for  the 
Ande r son- Da r 1 i ng  goodness-o f - f i t  test  (26).  The  computa¬ 
tional  form  used  in  this  case  would  be 

n2  •  -  %l>2- 


to 


(8) 


Problem  Statement 

Few  goodness-of- f i t  tests  are  available  to  perform  on 
sample  data  when  the  parameters  of  the  distribution  are  not 
Known.  As  mentioned  earlier,  Lilliefors  developed  a  test 
for  the  unspecified  exponential  and  normal.  Also,  bush 
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generated  a  set  of  critical  values  for  the  Weibull  with 
unspecified  scale  and  location  parameters  (3).  There  still 
exists  the  need  to  develop  a  valid  goodness-of-f it  test  for 
the  gamma  density  function  when  the  scale  and  location 
parameters  are  unknown. 

The  purpose  of  this  research  is  to  develop  a  goodness- 
of-fit  test  for  the  3  parameter  gamma  when  the  scale  and 
location  parameters  must  be  estimated  from  the  sample  data. 
This  involves  generating  a  table  of  critical  values  based  on 
the  sample  size  and  the  shape  parameter.  The  accuracy  of 
the  critical  values  must  be  sufficient  enough  so  that  data, 
sampled  from  other  populations  are  rejected. 

Objectives 

This  thesis  has  the  following  objectives: 

1)  To  generate  and  document  the  Anderson-Darling,  Crarner- 
von  rlisos,  and  the  Kolmogorov-Rmirnov  rejection  tables 
Lor  the  three  parameter  gamma  distribution  where  the 
scale  and  location  parameters  are  unknown. 

2)  To  conduct  a  power  comparison  between  the  Anderson- 
Darling,  Cramer-von  Mises,  and  Kol mogorov-Sm i rnov 
goodness-of-f it  tests. 

3)  To  investigate  the  possibilities  of  a  functional 
relationship  between  the  shape  parameter  and  the 
critical  values  in  objective  one. 


II.  The  Gamma  Di_s t r  i  but  i  on 
The  Cjmna  bent, i  ty  l‘'unc t ion 

The  gamma  density  function  is  useful  in  reliability  and 
Maintainability  theory.  It  also  has  applications  in  the 
natural  sciences.  If  the  random  variable  x  is  gamma  distri¬ 
buted  tnun  the  probability  density  function  taxes  the  form: 

f(x)  =  exp  ,  (9) 

r(K)0k  9 

9,  K  >  0;  x>  c  >  0. 

Ynore  are  three  parameters  which  specify  the  gamma.  0  is 
thu  scale  parameter;  K  is  the  shape  parameter,  and  c,  the 
location  parameter. 

Yhc  complexity  and  versatility  of  the  gamma  distribu¬ 
tion  can  be  observed  by  examining  the  graphs  of  the  distri¬ 
bution  for  various  values  of  the  shape  and  scale  parameters, 
figure  la  shows  graphs  of  the  standardized  gamma  (i.e.r  c=0, 
=  1)  tor  .  i>,  2 , 3 , 4  ,  and  5  (7:370).  When  K  =  l,  the  gamma  is 
Lite  exponential.  Also,  it  is  interesting  to  note  that  when 
is  loss  than  one,  the  gamma  closely  resembles  the  exponen¬ 
tial  distribution. 

/•lsu,  figure  lb  shows  that  for  large  K  (K=50),  the 
gam.  a  resembles  the  normal  distribution,  figure  2  illus- 
l  rates  tnc  influence  of  on  the  graph  of  tnc  jamma.  Here 
wo  set  c  =  uivi  \  -2  ana  sketen  the  graphs  for  0=1/3,  1  / 2 ,  1 


1  i  cat  ion  o  t  the  Gamine  Densi  ty  function 

The  gamma  distribution  is  oiten  used  in  reliability 
theory  to  represent  the  distribution  of  the  time  between 
Lai  lures  of  a  system.  Assume,  for  example,  that  a  system  is 
made  up  of  r  components,  all  of  which  must  fail  for  the 
system  to  fail.  furthermore,  assume  that  the  time  to  fail¬ 
ure  Xj  of  eacli  component  is  independent  and  exponentially 
distributed.  Then  the  time  to  system  failure  Y  =  +  X2  + 

...  +  Xr  is  gamma  distributed  (7:369). 

In  queueing  theory,  the  random  variable  T  follows  a 
Jumna  distribution,  where  T  =  Xi  +  Xp  +  •••  X  is  the  total 
time  to  service  K  customers  assuming  that  the  time  of 
service  of  each  customer  is  independent  and  exponentially 
ci  i  s t  r  i  Luted  . 

The  two  cases  described  can  be  modeled  as  a  Special 
case  of  the  .jamma  known  as  the  Erlang  distribution  and 

expressed  as: 

K 

f(t)  =  1_LL_  •  tK_1  e“K  fc,  t>0.  ( 1C ) 

( K-l )  !  " 

A  random  variable  having  a  gamma  distribution  has  also 
Leen  used  to  represent  or  measure  the  occurrence  of  physical 
phenomena.  lor  example,  black  and  Krubein  (1955) 
demonstrated  that  the  mean  value  x  of  radioactivity  (alpha 
particles  per  minute)  within  a  sample  of  Pennsylvania  shale 
followed  a  gamma  distribution  (7:37d). 
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Methodology 

This  chapter  presents  the  Monte  Carlo  simulation  proce¬ 
dure  used  in  generating  the  critical  value  tables  for  the 
modified  :<ol  nogorov-Sm  i  rnov,  Cramer-von  Mises,  and  Anderson- 
darling  tu'-sl:..  The  procedure  is  outlined  using  the  flow 
chart  in  figure  3.  Secondly,  an  outline  oi  the  power  com- 
parison  among  the  three  goodness-o f - f i t  tests  is  given. 
Thirdly,  a  discussion  of  the  analysis  of  the  functional 
relationship  between  the  shape  parameter  and  critical  values 
is  presented  for  each  of  the  test  statistics. 

Monte  Carlo  Simulation  Procedure 

The  following  procedure  is  used  to  generate  the  criti¬ 
cal  value  tables  for  the  modified  goodness-of-f i t  tests.  As 
mentioned  earlier  Figure  3  presents  these  steps  in  flow 
chart  format. 

1)  For  a  fixed  sample  size  n  and  fixed  shape  parameter  ;• , 
n  standard  random  gamma  deviates  are  generated  using  a 
computer  subroutine.  The  standard  gamma  deviates  are 
converted  to  random  deviates  with  location  parameter  C 
=  lb  and  scale  parameter  9  =  1. 

2)  The  n  random  deviates  are  ordered,  x(i)»  x(2)»  ••• 

X(n)  * 

3)  The  ordered  random  deviates  are  used  to 
estimate  the  maximum  likelihood  scale  and 
location  parameters. 
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4 )  The  estimated  scale  and  location  parameters  and  fixed 
shape  parameter  are  used  to  determine  the  hypothesized 


distribution  function  I'(x). 

9)  The  test  statistic  is  calculated  using  equations  three, 
six,  and  eight,  for  the  modified  Do  1 mogo r ov-Gm i r no v , 
Anderson-Darling,  and  Cramer-von  mises  tests  respec¬ 
tively. 

G)  Steps  one  through  five  are  repeated  5U0U  times.- 
7)  The  value  of  each  statistic  are  ordered  in  ascending 
order  and  the  80th,  89th,  90th,  95th  and  99th  percen¬ 
tiles  are  used  as  the  critical  values  of  the  test, 
donor  at  ion  o  t  the  Three  Parameter  Gamma  Deviates 

lor  the  gamma  distribution  function,  there  is  no  closed 
form  tor  which  we  could  obtain  an  inverse;  however  algo¬ 
rithms  are  available  which  can  be  used  to  generate  random 
gamma  deviates.  The  ImSL  subroutine  GGmAR  is  used  to  gener¬ 
ate  standard  gamma  deviates  in  this  thesis.  These  standard 
deviates  are  converted  to  deviates  having  location  C  =  10 
and  scale  0=1.  This  is  done  by  using  tne  transformation 

z  ■  •  •  x  +  C  (11) 

where  x  represents  a  standard  random  deviate.  This  trans- 
iori,  at  ion  is  made  to  avoid  a  problem  with  the  parameter 
estimating  routine,  further  discussion  on  this  matter  is 
presented  in  tne  following  section. 

i-md  iiiui:i  L  i  kel  i  hood  Dsti  mates  for  the  Gamma  Para  meters 

the  procedure  used  to  calculate  the  maximum  likelihood 
estimates  Lor  gamma  parameters  was  developed  by  liar  ter  and 


Moore  (12).  Their  analysis  involves  the  derivation  of  the 
maximum  likelihood  estimators  and  includes  an  iterative 


method  for  solving  the  simultaneous  equations.  To  derive 
the  maximum  likelihood  equations  we  begin  with  the  gamma 
density  1  unction  with  location  parameter  C>0,  scale  param¬ 
eter  9  ,  nriu  shat-e  parameter  K: 

f(x,c,#,K)  =  1  exp [-*j£]f  (12) 

.  9  ,  K>c  x>C>(). 

The  likelihood  function  of  the  order  statistics  x^,  x2,...xn 

o  l  a  s a m i  1  e  of  size  n  is 


L  =  ( _ A _ )n 

T(k)fl;  , 


)n  £  [it-p-l  exp  [-£i£d. 

i»l  *  i=l  * 


Ive  wish  to  find  the  values  of  9,  Kr  and  C  which  maximize  L. 
This  is  done  by  taking  the  natural  logarithm  of  L,  and 
setting  the  partial  derivatives  with  respect  to  the  three 
parameters  equal  to  zero  and  solving  the  three  simultaneous 
equations.  The  partial  derivatives  are  shown  here: 

dlnL  nk  Ji,  xi-c 

-sr  =  t  ♦  &  > 


-jr  -  -nln* +  £  ln0ti-c)-n  -%kKJ-  Tf 


rm 


dint 


=  (l-K)  £  (x.-c)-1  +  ~ 

i  =  l  • 


It  snould  be  noted  that  because  of  a  limitation  of  the 


barter  and  Moore  subroutine,  it  is  not  possible  to  estimate 
tut.  parameters  when  tne  gamma  deviates  are  generated  with 
location  C  =  C.  In  the  event  that  gamma  deviates  are 


t? 


cjenerated  with  C  =  t),  it  is  possible  to  obtain  negative 
estimates  for  the  parameters.  The  subroutine  maps  these 
negative  estimates  onto  zero.  Thus,  C  and  •  will  not  retain 
the  invariant  property  which  is  needed  for  these  tests  to  be 
valid.  In  addition,  the  iterative  technique  used  in  the 
i.arter  and  hoore  subroutine  does  not  work  for  the  special 
case  when  the  shape  parameter  is  set  to  one.  Because  the 
gamma  is  an  exponenti  n  1  c*  i  l.  t  r  i  i.>  ution  when  the  shape  param¬ 
eter  is  e-1ual  to  one,  we  can  use  the  maximum  likelihood 
estimators  of  the  location  and  scale  parameters  for  the 
exponential.  Therefore,  setting  k  equal  to  one  and  solving 
equations  13  and  lb  we  obtain 


C  =  x 


(1) 


ana 


1 1 

•  -  -1-  E 


n  ; - 


X:  -  X 


i  =  1 


(1) 


(17) 


(18) 


'  'c  r  i  v  i  ng  t_hc  llypothesi  zed  Distribution  1'unc  t  ion  !■'  ( x ) 

The  maximum  likelihood  estimates  Cor  the  location  and 

a 

scale  parameters  and  the  fixed  shape  parameter  determine  I- 
(x).  Obtaining  a  numerical  value  for  F(x)  requires  an 
integral  calculation;  this  calculation  was  done  using  the 
li'iSL  subroutine  iViDCAM  (13).  m  DC  AM  calculates  the  probabil¬ 


ity  that  a  random  variable  x  from  a  standard  gamma  distribu¬ 
tion  (i.e.,  C  -  i,  and  8=  1)  in  less  than  or  equal  to  x.  To 
t ranstor m  t:ao  deviates,  y^,  from  a  -jeneralized  gamma  distri¬ 
bution  into  standard  gamma  deviates  we  use 

a 

Vi-C 

Xi  =  -J— 


(lfJ) 


ir> 


The  details  oi  this  transformation  are  provided  in  Cortes 
(5: 17)  . 

Pool  r  Compel  r i son 

The  powers  of  the  modified  Kolinogcr  ov-Sm  i  r  nov ,  Cramer- 
von  oiises,  and  Anderson-Darling  tests  are  compared  for  ten 
alternative  distr ioutions.  Samples  of  sizes  equal  to  five, 
15,  and  2  5  are  drawn  from  the  following  selected 
d i s t  r i hot  ions: 


1) 

Gamma,  shape  equals 

1.5 

2) 

Gamma,  shape  equals 

2.5 

3) 

Gamma,  shape  equals 

4.1) 

4) 

Voibull,  shape 

equals  2. 

5) 

Ueioull,  shape 

equal 

s  3. 

f  ) 

Mori -a  1  ( 1  (! ,  1 ) 

7) 

beta  (p  =  1,  q 

-  2) 

P) 

i-,eta  (p  =  2,  q 

=  2) 

'■) 

Lognormal  (p  * 

l,  w 

=  0) 

Id) 

Lognormal  (p  * 

2,  w 

*  0) 

These  distributions  are  tested  according  to: 

H .. :  The  sample  variates  follow  a  gamma  distribution 

having  shape  parameter  K. 

The  sample  variates  follow  some  other  distribution. 
The  power  investigation  was  conducted  under  two  null  hypoth¬ 
eses,  one  Lor  shape  parameter,  =  1.5,  the  other  for  shape 
parameter,  K  =  4.'1.  The  rondo:.:  deviates  tor  the  above 

alternative  distributions  were  generated  using  IMSL 


subroutines 


The  lognormal  density  function  is  written  as 


t(;<)  = _ -1 _ exp  (-1  [liLJLZtfl  2  )  ,  :<>0  (  20 ) 

xpVzT  2  P 

-  .1  otherwise. 

and  is  ill  us  Lr<«  tec.  in  figures  4a,  and  4b  lor  the  parameters 
Pune  fa)  given  above. 

'me  nchiull  density  function  is 

•(>:)  =  liLl-d---1-  exp  t-iJiZ£i]K,  ;:,*>(!,  C^X  (21) 

flK  0 

=  ii  otherwise. 

and  is  stiown  in  figures  5a  and  5b. 

'i  he  beta  density  function  is  expressed  as 

f(x)  =  .-rO.^1-  x^d-x)1!’1,  t:<x<l  (22) 

=  (i ,  otherwise. 

nn<.  its  graph  for  the  two  cases  of  interest,  is  presented  in 
1  inures  Go  arid  6b. 

•  or  each  of  the  sample  s  i;;es  mentioned  above,  five 
thousand  so.ii-le  sets  were  generated  tor  the  alternative 
uistr i out  ions.  Vhe  location  and  scale  parameters  are  calcu¬ 
lated  under  the  null  hy(>othesis  and  the  three  test  statis- 
2  2 

tics,  I'-.:,  A  ,  and  W  arc  evaluated.  The  value  of  these 
statistics  ate  compared  to  the  critical  values  derived  in 
this  t.iesm.  il  the  value  of  the  statistic  is  greater  than 
the  critical  value,  the  null  hypothesis  is  rejected.  The 


■ 


1 U 


total  number  of  rejections  are  counted.  The  power  is  the 
total  number  of  rejections  divided  by  the  number  of  trials, 

f.bMO. 

Determining  the  Critical  Values 

t.y  repeating  steps  one  through  five  5U0U  tines  as  shown 
in  the  flow  chart  in  figure  7,  5f!0!  values  for  K-S,  A  ^ ,  and 
Vt  are  calculated.  These  critical  values  are  ordered  and 
the  Oi/th,  85th,  90th,  95th,  and  y9th  percentile  are  used  as 
the  critical  values  of  tne  tests. 


Computer 


Protj  rams 


The  computer  programs  used  in  this  thesis  are  presented 


in  Appendix  G. 
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I V .  Use  of  Ta b  1  e  s 

In  this  chapter  a  set  of  steps  is  given  which  is  used 
to  pert'ori;i  a  goodness-o i- f i t  test  by  applying  any  of  the 
three  tests  developed  in  this  thesis.  Also,  an  example 
illustrating  the  Anderson-Darling  test  is  presented. 

The  folowing  steps  are  used  to  perform  a  goodness-of- 
L i t  test: 

1)  Determine  tne  shajpe  parameter,  K,  and  the  desired  level 
uL  significance  a  . 

2)  i  roci  the  data  to  be  tested,  calculate  the  maximum 
likelihood  estimators  for  the  location  and  scale 

t  )ci  T  I  kU  Let  o  • 

3)  iron  the  appropriate  table,  select  the  critical  value, 

dft  a  corresponding  to  <*  ,  the  sample  size  n,  and  shape 

parameter  K. 

4)  Using  Lite  maximum  likelihood  estimators,  determine  the 
estimated  hypothesized  distribution,  and  use  equation 
three,  six,  or  eight  to  calculate  the  Kolmogorov- 
duirnov,  Anderson-Darling,  or  Cranor-von  Mises  test 
statistic  respectively. 

h)  II  tin:  value  obtained  in  step  four  is  greater  than  the 

critical  value  found  in  step  three,  then  reject  the 

hypothesized  distribution.  11  it  is  smaller  than  the 
critical  value,  then  the  hypothesized  distribution  can 
not  ./v-  rejected. 


I 
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LI  X iM in i  >  It' 


The  time  between  failures  of  a  particular  subsystem  of 
a  radar  system  is  believed  to  be  distributed  according  to  a 
gamma  distribution  with  shape  parameter  equal  to  3.0.  A 
test  engineer  recorded  the  following  times  between  failures 
of  that  subsystem:  11.1,  10.6,  10.4,  1 3 .  (> ,  11.3,  10.5, 

1C. 6,  10.9,  10.6,  10.8  days. 

/•.  modi  Lied  Anderson -Darling  test  at  a  .0/5  level  of 
sign  i  l  icance  is  performed  using  the  critical  values  in  this 
tnesis.  The  problem  can  be  stated  as  a  test  of  hypothesis, 
that  is, 

il^ :  The  distribution  is  gamma  (shape  =  3.0). 

ii A :  The  sample  comes  from  another  distribution. 

first,  the  level  of  signi  f  icance  o  ,  and  shape  parameter  K, 

have  been  determined  to  be  .05  and  3.0  respectively. 

second,  the  maximum  likelihood  estimators  calculated  using 

tne  darter  and  moore  subroutine  are  d  =  9.319  and  •  =  .026. 

Next,  the.  critical  value  from  Table  XXVI  is  .8415.  The 

hypothesized  distribution  is  completely  determined  by  the 

fixed  shape  parameter  and  the  estimated  location  and  scale 

parameters;  these  values  are  presented  in  Table  I.  The 

.  .  .  .  0 

value  ol:  the  And e r son- Da r 1 l ng  statistic  is  A  =  1.734  2. 
Since  1.7342  is  greater  than  .0415,  the  null  hypothesis  is 
r ejected.  Therefore,  the  conclusion  is  that  the  sample  of 
time  between  failures  comes  from  some  otner  distr.U.  tion. 


V .  Discussion  oi  the  tfesu 1 ts 

This  chapter  presents  the  results  obtained  with  respect  to 
the  objectives  stated  in  chapter  1.  These  objectives  were  to 
develop  modified  Kolmogorov-Smirnov,  Anderson-Darling,  and 
Cr iiuer -von  Hisos  tests  tor  the  gamma  ana  compare  their 
powers.  Also  included  was  an  investigation  of  the  relation¬ 
ship  between  tne  critical  values  of  each  test  and  shape 
parameters.  Included .with  these  results  is  a  report  on  the 
validation  of  the  computer  programs  used  in  this  thesis. 

Presentation  of  tne  Kolmogorov-Smirnov,  Cramer-von  Mises, 
and  Aimer son-ma r 1 i ng  Tables  of  Cr i t i ca 1  Values 

The  tables  of  critical  values  for  the  modified  K-S,  A  , 
ana  ^  tests  are  presented  in  Appendices  A,  B,  and  C, 
r  espuc L i vel y . 

o 

‘..hen  t. he  snape  parameter  is  fixed,  both  the  K-C  and  A 
critical  values  are  decreasing  as  the  sample  size  increases. 
The  rate  ot  decrease  is  smaller  as  n  increases;  this  is  an 
ir. dication  that  the  critical  values  appear  to  be  converging 
for  large  sample  sizes.  The  Cramer-von  mises  critical 
values,  on  the  other  hand,  are  increasing  with  respect  to 
the  sample  size.  Again,  the  rate  of  increase  is  smaller  for 
larger  values  of  n,  indicating  that  the  critical  values  are 
converging  for  large  sample  sizes. 

It  should  be  noted  that  because  the  critical  values  are 
derived  through  Monte  Carlo  simulations,  the  values  are  not 
error  1  rue  and  tnat  the  ,  amount  of  error  decrease's  as  the 
numi.or  of  trials  increases  (25).  The  bhtib  repetitions  used 
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in  this  thesis  was  a  practical  compromise  based  on  computer 
time  required  and  accuracy  desired. 

Computer  Programs  Validation 

The  computer  programs  are  verified  by  generating  crit- 

* 

ical  values  for  the  exponential  distribution  with  unknown 
mean.  This  is  done  by  generating  gamma  deviates  with  shape 
parameter  equal  to  one,  fixing  the  location  parameter  at 
some  arbitrary  value,  and  estimating  only  the  scale 
parameter 

The  critical  values  are  calculated  for  sample  sizes  n  = 
5,  lb,  20,  and  30.  The  critical  values  from  the  Anderson- 
darling  arid  Craraer-von  Mises  statistics  are  modified  using 
expressions  (23)  and  (24)  derived  by  Stephens  (27): 

A2  (1  +  -  -|)  (23) 

n  n2 

a2  (1  +  ^iil)  .  (24) 

n 

The  commuted  critical  values  are  compared  to  those 
calculated  by  Stephens  (27)  for  significance  levels  .15, 

Vue  critical  values  calculated  for  the  Ko 1 mogo r ov-Sm i r nov 
statistic  for  the  exponential  are  compared  ditectly  to  those 
derived  by  Lilliefors  (20).  The  critical  values  which  are 
derived  from  the  programs  in  this  thesis  are  presented  in 
Table  IV  and  can  be  compared  to  Lilliefors  results  in  Table 
V  . 

The  critical  values  for  the  Kol mogo rov-Sm i rnov  statis¬ 
tic  compared  very  well  to  the  Lilliefors  values.  The 
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TABLE  II 

.  ■) 

Cramer- von  Mi sea  W 


w2  (1  +  ,16/n) 


.  140 
.174 
.  221 
.335 


.1523 
.  17C 
.  221 
.325 


Stephen ' a 
Critical 
Values 


.  149 
.171 

.  224 
.337 


TABLE  III 

Anderson-Dar 1 in j  A2 


A2  (1  +  1.5/n  -  5/n") 

Stephen '  s 

*  X 

n=b 

- - 

3 

II 

•-> 

n  =  2C 

n  =  30 

c  r 1 1 1 ca i 
Values 

2  5 

.941 

.992 

.953 

.901 

.922 

5  2 

1 . 955 

1.129 

1.097 

1.098 

1.070 

5  5 

1.357 

1.351 

1.37C 

1.372 

1.341 

99 

2  .  >’  /  9 

2.219 

2.077 

1.974 

1.957 

greatest  deviation  occurs  for  n  =  5  at  significance  level 
.01.  The  Cramer-von  Mises  values  generated  are  very  close 
to  Stephens  values  with  the  greatest  deviation  being  3.Gfc 
for  n  =  3U  and  significance  level  at  .01.  There  was  also  a 
good  match  between  the  Ander son-Da  r  1  i  ni  values,  with  most 
deviations  between  3't  and  4%;  however  the  greatest  deviation 
is  13. -i0  for  n  =  10  and  significance  level  . i'l. 
bower  Irivestigatio  n 

n  [lower  co m pa r  i  so n  is  made  between  the  Kolmogorov- 
Smirnov,  Anderson- Darling,  and  Cramer-von  Mises  goodness-of- 
L'it  tests  developed  in  this  thesis.  The  gamma  distribution 
with  shape  parameters  equal  to  1.5  and  4.(5  were  both  used 
against  the  alternative  distributions  listed  in  chapter  III. 
Sample  sizes  five,  15,  and  25  were  used  in  the  power  studies 
at  both  an  <*-level  of  .05  and  .01. 

Tables  VI  through  XIII  show  the  results  of  the  power 
comparisons  tor  a -levels  .05  and  .01.  When  tire  null  hypoth¬ 
esis  is  true,  the  power  meets  the  claimed  level  of  signifi¬ 
cance  to  the  second  decimal  place  in  most  cases.  For  all 
test.,  the  power  is  low  for  sample  sizes  equal  to  five.  In 
tact,  in  most  cases  for  n  =  5  the  power  is  nearly  equal  to 
the  significance  of  the  test,  indicating  that  the  goodness- 
of-fit  test  has  no  practical  use  for  very  small  sample 

*  *  i  Z  G  iS  • 

In  ner.rly  ail  cases,  the  powers  of  the  Cramer-von  Mises 
and/or  anaer son-Dar 1 ing  are  greater  than  Kol mogorov-Smi rnov 
tests.  based  on  this  study,  the  latter  test  would  not  be 


istr ibutior 


for 


bower  Teat  tor  the  Gamma  Distribution 
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used  as  long  as  the  first  two  are  available.  An  examination 
of  Tables  VI  through  XIII  reveals  that  for  both  levels  of 
significance,  the  Cramer-von  Mises  test  is  more  powerful 
than  the  Anderson-Darling  test.  The  only  exception  to  the 
previous  statement  occurs  when  the  alternative  distribution 
is  lognormal. 

The  following  observations  are  made  concerning  the 
power  of  both  the  Cramer-von  Mises  and  Anderson-Darling 
tests  when  sample  sixes  are  either  15  or  25: 

1)  It  the  null  hypothesis  is  a  gamma  with  shape  equal 
to  1.5,  too  power  is  high  against  all  alternative  distribu¬ 
tions  except  tor  another  gamma.  The  tests  are  especially 
high  against  the  lognormal  distribution. 

2)  When  the  null  hypothesis  is  the  gamma  with  shape 
equal  to  four,  the  power  is  high  against  the  lognormal  and 
normal  distributions.  The  power  is  not  qu;te  so  high 
against  the  beta.  Against  the  other  gammas  and  the  Weibull 
with  shpae  equal  to  two,  the  power  is  low,  even  when  the 
sample  size  is  25. 

Kola  L  i  onsh  i  [■  me  tween  Critical  Values  and  Shape  Parameters 

An  investigation  of  the  relationship  between  the  shape 
parameter  and  the  Kolmogorov-Smirnov,  Anderson-Darling,  and 
Cramer-von  wises  critical  values  is  summarized  in  this  sec¬ 
tion.  For  each  test  statistic,  the  shape  parauieter  versus 
critical  values  are  plotted.  These  graphs  arc*  presented  in 
Appendix  D  Lor  the  K-S,  Appendix  E  for  A2,  and  Appendix  F 
for  the  2  critical  values.  The  graphs  of  all  the  test 
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statistics  appear  to  exhibit  a  common  and  consistent  behav¬ 
ior.  This  consistent  behavior  observed  from  the  graphical 
representations  can  be  summed  up  as  follows: 

1)  tor  all  test  statistics,  the  relationship  of  crit¬ 
ical  values  as  a  function  of  shape  is  always  decreasing  for 
shape  greater  than  one;  this  decrease  appears  to  be  an 
inverse  relationship. 

?)  The  graphs  of  the  Anderson-darling  critical  values 
always  show  an  increase  as  the  shape  increases  from  .5  to 


1 .  0 . 


3)  The  holmogorov-Smirnov  and  Cramer-von  Mises  crit¬ 
ical  values  increase  or  decrease,  as  the  shape  varies  from 
.b  to  1.0,  aepending  on  the  sample  size. 

A  regression  analysis  is  performed  to  determine  the 
functional  relationship  between  the  critical  values  and 
shape  parameters,  as  suggested  by  the  graphs.  The  study 
includes  values  for  the  shape  between  1.5  and  4.0.  Shape 
parameters  less  than  1.5  are  not  considered  in  the  regres¬ 
sion  analysis  because  a  different  estimating  technique  was 
used  to  calculate  the  critical  values  for  the  gamma  when 
shu is  e  4 u a  1  to  one.  In  addition,  more  information  is 
needed  about  the  behavior  of  the  function  between  .5  and  1.5 
for  ull  test  statistics. 

The  expression  which  best  represents  the  relat ionsh i p 
Lor  all  test  statistics  is 

C  a  a^,  +  a  j  ( 1/K'  )  . 

k‘  ,  the  number  whicn  measured  the  amount  of  variation 
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in  ail  cases.  This  expression  can  be  used  to  find  the 
critical  value  corresponding  to  shape  parameters  between  1.5 
and  4.U,  not  found  in  the  tat>les  in  Appendices  A,  li,  and  C. 
Therefore,  a  test  of  hypothesis  can  be  performed  for  the 
null  hypothesis  being,  lor  example,  a  gamma  distribution 
having  shape  equal  to  2.7  5. 

The  value  of  K  2 ,  and  coefficients  a(,  and  a^  are  re¬ 
corded  in  Tables  XIV,  XV,  and  XVI,  for  the  K-S,  A2,  and  Vn,  2 
critical  values  respect ively. 
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VI.  Conclusions  and  Hecomncndations 

Concl usions 

based  on  results  obtained  in  this  thesis,  the  following 
conclusions  are  noted: 

1)  The  Kolmogorov-Smirnov,  Anderson-Darling  and 
Cramer-von  Mines  critical  values  for  the  three-parameter 
gamma  are  valid.  The  power  study  revealed  that  when  the 
null  hypothesis  is  true  all  three  tests  achieve  the  claimed 
level  of  significance. 

2)  The  power  comparison  study  based  on  the  ten  alter¬ 
native  distributions  listed  in  Chapter  3  shows  that  in 
general  the  powers  in  decreasing  order  are  W2,  A2,  and  K-S. 
The  A2,  however  is  more  powerful  against  the  lognormal 
distribution.  All  throe  tests  demonstrated  low  power  for 
sample  sizes  ogual  to  five,  indicating  a  goodness-of-f i t 
test  involving  a  sample  size  of  five  using  the  tabled 
critical  values  would  not  be  practical. 

i  ’ccoi  i.  ,e nd a  t  ions 

The  following  recommendations  are  suggested  for  further 
i nves t igat ion : 

1)  Develop  a  more  efficient  technique  to  calculate 
the  maximum  likelihood  estimators  for  the  parameters  of  a 
gamma  distribution. 

2)  Investigate  a  functional  relationship  between  the 
sample  size  and  critical  values  so  that  goodness-o f-f i t 
tents  can  be  done  for  sample  sizes  other  than  those 
presented  iri  this  thesis. 
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3)  t’xtend  the  goodness-of-fit  test  to  include 
parameters  between  zero  and  one. 

4)  Examine  the  feasibility  of  developing  a  yoodness- 
of-fit  test  for  distributions  whose  parameters  are  unknown, 
based  on  the  characteristic  function. 
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TABLE  XIII 


Kolmogorov-Smi rnov 
Shape  Parameter  =  .5 


S  an  p  1  e 

Level 

of  Significance 

n 

.  20 

.15 

.10 

.05 

.01 

t; 

V 

.3528 

.  3681 

.3913 

.  4330 

.5217 

10 

.2537 

.2798 

.3010 

.3331 

.3862 

15 

.  2249 

'  .  2365 

.2514 

.  2804 

.3317 

20 

.1967 

.2077 

.2225 

.2457 

.2941 

25 

.  17  59 

.1870 

.2006 

.2227 

.2623 

30 

.1635 

.1736 

.1844 

.2(139 

.2449 

TABLE  XIV 


Kolmo'jorov-Smi  rnov 
Shape  Parameter  =  1.0 


Sam pie 

Level 

of  Sicjni f icancc 

J  1  /<  L 

n 

.  20 

.15 

.10 

.05 

.01 

5 

.3701 

.  3848 

.  3994 

.  4433 

.5269 

10 

.2651 

.2788 

.2958 

.3267 

.3930 

15 

.2184 

.  2  306 

.  2451 

.  2709 

.  3203 

O  «• 

/_  tj 

.1923 

.2038 

.2177 

.239  2 

.2774 

25 

.  1693 

.1788 

.1910 

.2113 

.  2520 

30 

.1561 

.1650 

.17  55 

.1958 

.2297 

ri  pI  e 


TABLE  XVII 

Kolmogcrcv-Smi rnov 
lliape  Parameter  =  2 


Level  of  Significance 


5 


3198 

.  3  366 

.  3581 

.3898 

2370 

.2483 

.2634 

.2884 

1963 

.  2867 

.2187 

.2414 

1718 

.1803 

.1911 

.2102 

1552 

.1636 

.  1738 

.1890 

1417 

.1483 

.1584 

.1730 

.4399 
.3378 
.2814 
.2457 
.  2289 
.2023 


TABLE  XVIII 

Kolmogorov-Smirnov 
hape  parameter  =  3 


Sample 

Level 

ot  Significance 

s  i  ze 
n 

.20. 

.15 

.10 

.85 

.01 

1 

.3179 

.  3342 

.  3524 

.  3799 

.4  365 

18 

.2338 

.  2449 

.259  5 

.28  29 

.3282 

1‘. 

.  1939 

.  2041 

.2172 

.2352 

.2751 

20 

.  1710 

.1792 

.1898 

.2067 

.2395 

25 

.1520 

.  1600 

.  1695 

.  1858 

.2130 

3(. 

.1411 

.1478 

.1573 

.1711 

.2026 
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TAD  Li!  XIX 


Kolmoyorov-Smi r nov 
SLa^e  Parameter  -  3.5 


Level  of  Significance 


.20 

.15 

.10 

.05 

-01 

3161 

.3317 

.3508 

.  3769 

.4314 

2310 

.2421 

.2561 

.2812 

.3234 

1923 

.  2001 

.2127 

.  2326 

.  2724 

169  5 

.1775 

.1084 

.20  59 

.2403 

1511 

.1580 

.  1689 

.  1029 

.2127 

138  2 

.1452 

.1543 

.1604 

.199  2 

TAIiLE  XX 


Kolmogorov-Smi rnov 
Shape  Parameter  =  4.',' 


Level  of  Significance 


.  20 

.15 

.10 

.05 

.01 

,  3131 

.  3289 

.3471 

.  3731 

.4266 

.2300 

.2416 

.2567 

.2792 

.3267 

,  1904 

.  1992 

.2116 

.2312 

.  2666 

.  1667 

.1740 

.10  58 

.200  2 

.2312 

.1505 

.1570 

.1666 

.  1023 

.  2137 

.1381 

.1450 

.1544 

.1690 

.1970 
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TABLE  XXI 


Ander son- Dari  inn 


Shape 

Parameter 

=  .5 

Level 

of  Signi ficance 

.  20 

.15 

.05 

.01 

1.0724 

1.1679 

1.3135 

1.6238 

2.5899 

.5224 

•  1.0295 

1.1919 

1.4375 

2.2323 

.8797 

.9824 

1.1258 

1.4083 

2.2116 

CO 

• 

.9759 

1.1352 

1.4831 

2.1989 

.8367 

.9472 

1.0963 

1.3910 

2.0626 

.8332 

.94  28 

1.1107 

1.4316 

2.3048 

TABLE  XXII 
Anderson-oarl  in<j 


shape 

Parameter 

=  1  .  fl 

Level 

of  Significance 

.  20 

.15 

.10 

.05 

.01 

2.0880 

2.3079 

2.5926 

3.0710 

4.3921 

1.5085 

1 . 6604 

1.8807 

2.2980 

3.4906 

1.2827 

1.4396 

1.6493 

2.0360 

2.9562 

1.21 2 4 

1.3510 

1.5480 

1 .915  38 

2.7789 

1  .'1831 

1.2147 

1 . 3854 

1.7118 

2.5678 

1.0558 

1.1750 

1.3756 

1.6694 

2.4482 

TABLL  XXIII 


Anderson-Darling 
Shape  lJaraiiieter  =  1.5 


Level  of  Significance 


mm 

.15 

■BOH 

.05 

.01 

GO  18 

.7  577 

.8451 

1 . 0603 

1.430C 

7104  , 

.7032 

.9017 

1.0520 

1.5597 

7063 

.794B 

.9007 

1.1067 

1.5828 

714  5 

.7903 

.9120 

1.0523 

1.5719 

,7278 

.8053 

.9170 

1.1332 

1.5569 

,7  7.96 

.8105 

.9266 

1.1189 

1.6133 

TABLL  XXIV 


Ander son- Dari i ng 
Shape  Parameter  -  2.0 


r 

Level 

of  Significance 

.20 

HREH 

.10 

.05 

.01 

.6106 

.6666 

.7417 

.8673 

1.1757 

.6  265 

.6949 

.7879 

.9240 

1.2260 

.6  330 

.6927 

CO 

• 

.9405 

1.3451 

.64  22 

.7101 

.8051 

.'/.V)  2 

1 . 3437 

|  .6429 

.7002 

.7984 

.  9  670 

1.3956 

.'>481 

.7167 

.8213 

.9976 

1.4640 

TABLE  XXV 

Anderson-Dar 1 ing 
hape  Parameter  = 


Lam ole 

S  i  ze 

Level 

of  significance 

n 

.20 

.15 

.10 

.05 

.01 

5 

.  5700 

.6261 

.6964 

.8324 

1.0970 

.  soo  i 

.6504 

.74  09 

.8740 

1.2080 

1  15 

.5015 

.6556 

.7330 

.8693 

1.2046 

2'.) 

.6109 

.6699 

.7504 

.9G29 

1.2276 

25 

.5991 

.6590 

.7462 

.9023 

1.2320 

30 

1 

.5011 

.6573 

.7394 

.9025 

1.2614 

TABLE  XXVI 


Anderson- Dari ing 
shape  Parameter  = 


Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

.5549 

.596  2 

.6730 

.7013 

1.0634 

.5607 

.6276 

.70  56 

.  8  415 

1.1702 

.  5757 

.6376 

.7251 

,0615 

1.1651 

.5749 

.6360 

.116' 

.0755 

1.1499 

.5700 

.6364 

.7225 

.0  808 

1.1714 

.5852 

.6490 

.7355 

.0028 

1.2601 

.20 

.15 

.  In 

.05 

.01 

5293 

.  5P  35 

.6  464 

.7535 

.9715 

547  5 

.5996 

.6720 

.8031 

1.1020 

54  35 

.  6  0  0  6 

.6816 

.8154 

1.1036 

5  50  5 

.61)60 

.6915 

.8204 

1.1021 

5557 

.6124 

.6998 

.8425 

1.2094 

557  7 


G19  3 


69  7  O 


8  239 


l.ir.iii 


TAbLE  XXIX 


Cramer -von  Mises 
hape  parameter  = 


Level  of  S  i.jn  i  f  icance 


1217 

.1363 

.1576 

.  1928 

1327  • 

.1530 

.179  5 

.2259 

1379 

.1596 

.1063 

.  2340 

1425 

.1625 

.1910 

.2526 

1453 

.1654 

.1951 

.2514 

1440 

.1649 

.1986 

.2568 

.3200 
.3569 
.  3660 
.3743 
,3769 
4170 


TABLE  XXX 

Cramer-von  Mises 
Shape  Parameter  = 


Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

.1397 

.1561 

.1009 

.2314 

.  3571 

.1341 

.1511 

.1775 

.2210 

.3573 

.1310 

.  1490 

.  1730 

.2192 

.  3320 

.1316 

.1484 

.1748 

.2220 

.3327 

.  1299 

.1487 

.  1742 

.2250 

.3239 

.1319 

.1501 

.1757 

.2220 

.3415 
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table;  xxxi 


Cramer-von  Kises 
shape  Parameter  = 


-ample 

Size 

pi 


Level  ol  Significance 


.1741 
.1015 
.  1846 
.1917 
.1847 
.1870 


.20 

.15 

.10 

1092 

.1210 

.1390 

1132 

.1268 

.1471 

1138 

.  1287 

.1541 

1199 

.1353 

.1551 

1157 

.1307 

.  1503 

1171 

.1324 

.1520 

2  36  3 
2670 
2740 
,2902 
,2791 
,2872 


T/'BLE  XXXII 

Cramer-von  Mi. yes 
hape  Parameter  =  2. 


‘..ample 
s  i  /,  c* 

Level 

of  Significance 

n 

.  20 

.15 

.10 

.05 

.01 

5 

.  1007 

.1118 

.1261 

.  1481 

.2096 

1  v; 

.  1045 

.1160 

.1336 

CM 

VT 

• 

.2228 

15 

.  1840 

.1154 

.1325 

.  16  27 

.  2409 

20 

.1055 

.1169 

.  1355 

.  16  8  6 

.2404 

25 

.  1051 

.1174 

.1356 

.  1695 

.  2481 

30 

.  1077 

.1228 

.1415 

.1736 

.2590 

Tam  Lb  XXaIXI 


Cramer -von  i*iises 
Ua.je  I’arameter  = 


Level 

of  Significance 

.20 

.  Id 

.05 

.01 

.095-3 

.  1053 

.1137 

.1/532 

.  1901 

.0309  • 

.1033 

.125  3 

.1522 

.2194 

.  0933 

.1100 

.  1251 

.1531 

.2207 

.0991 

.1097 

.1256 

.1531 

.2207 

.09  79 

.1099 

.  1277 

.1547 

.  2209 

.0973 

.1091 

.1255 

.1543 

.2237 

TABLE  XXXIV 


Cramer-von  Miseii 
iiape  parameter  =  3 


;  Oiii'i  pi  e 

L  i  Z>: 

n 

Level 

of  Significance 

.20 

.15 

.10 

.05 

.01 

r 

> 

.0927 

.  1011 

.1145 

.  1  367 

.  1398 

10 

.  094  0 

.  1042 

.12.10 

.14  11 

.2052 

15 

.09  54 

.1075 

.1234 

.1501 

.  2085 

?  ' 

.994  2 

.1057 

.1214 

.1514 

.  206  0 

25 

.0935 

.  1050 

.1224 

.  1524 

.2137 

3-/, 

.0,957 

.1031 

.1247 

.1521 

.  228  0 

TAULL  XXXV 


Cramer -von  Mises 
hape  Parameter  =  3.5 


Level  of  Significance 


.20 

.15 

.  10 

.05 

.(51 

.0919 

.  1009 

.1127 

.1341 

.1827 

.0914  ■ 

.1006 

.1148 

.1396 

.2007 

.0925 

.1032 

.1172 

.  1410 

.  2002 

.09  39 

.1051 

.1203 

.1473 

.2112 

.0920 

.10  34 

.  1189 

.1447 

.2116 

.0926 

.1(524 

.1170 

.1416 

.2056 

TABLL  XXXVI 

Cramer-von  Mises 
Shape  Parameter  =  4.0 


Level  of  Significance 


.  20 

.15 

.10 

•  5 

.01 

.0839 

.0968 

.  1109 

.1295 

.1763 

.0906 

.10  0  2 

.1150 

.1391 

.1974 

.0915 

.1012 

.1152 

.  1396 

.  1967 

.0908 

.  1015 

.1158 

.1424 

.19  59 

.0916 

.  1  (5  2 1 

.1182 

.  1433 

.2110 

.09  28 

.1019 

.1178 

.14  30 

.2050 

SHAPE  VS  CRITICAL  VALUES 
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GRMMR 


1.00  2.00  3. 

SHAPE  PARAMETER 


SHAPE  VS  CRITICAL  VALUES 


L-'lcj  52.  l.liaj'C!  vs  A2  Critical  Values,  Level  =  .IV,  n  =*  20 

114 


~r  r  r 


SHAPE  VS  CRITICAL  VALUES 


3 


CRITICRL  VALUES 

.60  0.80  1.00  1.20  1.40  1.60  I.  BO 


SHAPE  VS  CRITICAL  VALUES 
GAMMA 

LEVEL-. 15  N=15 


SHAPE  VS  CRITICAL  VALUES 


i  IU  A L. .  vt.  A2  Critical  Values,  Level  =  .20,  n  =  2'j 


127 


SHAPE  VS  CRITICAL  VALUES 


0.20 


SHAPE  VS  CRITICAL  VALUES 


SHAPE  VS  CRITICAL  VALUES 


SHAPE  PARAMETER 


il,  74 .  v.,  y.  Critical  Values,  Level  =  .LS,  n  =  1U 


137 


CRITICAL  VALUES 

.14  0.16  0.18  0.20  0.22  0.24  0.26 


SHAPE  VS  CRITICAL  VALUES 
GAMMA 

LEVEL=. 05  N= 15 


138 


i'.OO 2 '.DO  3 . 00 

SHAPE  PARAMETER 


hape  vs  U  Critical  Values,  Level 


Critical  Values,  Level 


SHAPE  VS  CRITICAL  VALUES 


.00 


1.00  2.00  3.00 

SHAPE  PARAMETER 


CRT 


ICflL  VALUES 


l  IG  L7 .  r.nal  C  vs  l\2  Critical  Values,  Level  =  .15,  n  = 


151) 


SHAPE  VS  CRITICRL  VRLUES 


SHAPE  VS  CRITICAL  VALUES 


Level 


Critical  Valuer, 


CRITICAL  VALUES  *1CT 

JD.92  1.00  1.08  1.16  1.24  1.32  1.40 

cr  i  i  i  _i_  i  ±  _ i 


SHAPE  VS  CRITICAL  VALUES 
GAMMA 


SHAPE  PARAMETER 


i3.  Jha^o  Vo  a  Critical  Values,  Level  =  .20,  n  *  15 


15f> 


SHAPE  VS  CRITICAL  VALUES 


PROGRAM  AO 


c  . . . . . . . 

C  ‘THIS  PRGGRAM  GENERATES  THE  MODIFIED  A-D  STATISTICS  • 

C  *5000  REPS  • 

C  THE  TABLES  GENERATED  ARE  VALID  FOR  THE  GAMMA  OISTR IBUTI OH  • 

C  • N  =  SAMPLE  SIZE  =  5. <51,30  • 

C  • SSI  s  t  IF  SCALE  PARAMETER  (THETA)  IS  KNOWN  • 

C  • SSI  :  1  IF  THETA  IS  TO  BE  ESTIMATED  • 

C  •  SS2  s  8  IF  SHAPE  PARAMETER  (K>  IS  KNOWN  • 

C  • SS2  =  1  IF  K  IS  TO  BE  ESTIMATED  • 

C  •  SS3  =  0  IF  LOCATION  (C)  IS  KNOWN  • 

C  *SS 3  =1  IF  C  IS  TO  BE  ESTIMATED  • 

C  *C1  =  INITIAL  ESTIMATE  OF  C  (OR  KNOWN  VALUE)  • 

C  • T1  =  INITIAL  ESTIMATE  OF  THETA  (OR  KNOWN  VALUE)  • 

c  *ai  =  initial  estimate  of  alpha  (or  known  value)  • 

c  ••••••••*•••*•••••••.••••••••••*•••••••*«•••••••*••••* . •  •••• 

c 

COMMON/VALUE /Pit  00) 

COMMCN/R A Y/T (100) 

COMMON/MANA/N* SSI *SS2* SS3* M* Cl • T 1 «A1 «MR 
DOUBLE  PRECISION  OSEEOt T »CI » T I  * Al 

DIMENSION  FX(60)*AA(SOQO)(XX(S002)*VT(5002>(6(SO) 

integer  rep.pp 
OSCE 0=150  00*000 
MR-0 
NON£=0 
NZERO-O 
REP- 1>0  02 
N0S=REP“2 
MUM=REP-2 
C 

c 

C  ‘CALCULATES  5000  PLOTTING  POSITIONS  ON  INTERVAL  (0«I)  AT  (I-.5)/N 

c  ‘These  are  useo  for  interpolations  on  percentiles 

c 

VV(l)=0. 
rV(REP): I* 

00  AOS  L-2, REP-1 

VY(L)=((L-1)-.S)/N0S 
AOS  CONTINUE 

C  • ENO  OF  PLOTTING  POSITIONS  ROUTINE 
READ*,SSl,SS2,SS3«Cl,Tl«Al 
PRINT *,SS1,SS2 ,SS  3«C 1 , T 1 , A l 
PRINT* 

PRINT* 

PRINT* 

PRINT*(2X,A«FS*l)*,*SHAPE=*,Al 
PRINT  »  <2* , A1 •  ,*— — — 

PRINT* 

C  ‘SAMPLE  SIZE  ASSIGNED  HERE 

00  100  PP=1S»IS 
N=PP 
M=N 

C  ‘LOOP  FOR  MONTE  CARLO  SIMULATION  STARTS  HERE 

00  99  KK=1,SOOO 
CALL  GGAMR(0SEE0f Al.N.GtP) 


1G1 


00  719  I K= 1 • N 
.  PIIK):l. •P(IK»«10. 

719  CONTINUE 

CALL  VSR  T  A (P.N) 

00  3  11=1. N 
T ( I I >  =  PC III 
3  CONTINUE 

CALL  CANNA(CSJ«rSJ(ASU> 

C  • CALCULATE  S  ESTIHATEO  FIX*  FOR  EACH  SAMPLE  POINT 
00  333  L  =  1  .N 
U=CP(L»-CSJI/TSJ 
X1=ASJ 

CALL  HOGAHIM.Xl.PROB.lERl 

FXCL )=PROB 

IFtFXCLI.Ca.O.irHEN 

FX(L>=FX<L>*.0001 

N2£R0=NZER0*1 

END  IF 

IF<FX<L).£Q.1.)THEN 
FX«LI=FX(L>-.Q0Q1 
N0NE=N0NE*1 
END  IF 

333  CONTINUE 
UA0=0. 

XN-N 

00  300  I=1,N 
XI=1 

UAO=UAO*<2.*XI-l.}*iLOG<FK(I )>*L0G(1.-FX(N*1-I>n 
500  CONTINUE 

MAD=<-HA0/XN)-XN 
AA«KK)=UAO 
99  CONTINUE 

CALL  VSRTATAA.SOOOl 
00  400  L  =  l, REP-2 
XX(L*1 )= AACL I 
400  CONTINUE 

CALL  ENDPTTXX.YV .REP. NUN) 

C 

C  ‘PRINTS  PERCENTILES 

C 

PRINT((2X«AtI2)'«*F0R  N  =  *,PP 

PR  IN T*  (2  X  t  A)  »,  • - - - • 

PRINT* 

00  410  J  =  80  «  9  5  .5 
00  420  I  1  =  1. REP 
I=REP*1- 1 1 

IF«VT(II.LT.(J/100.0)ITHEN 
SL0PE=<T»(I‘1)-TT(I ))/(XX(I»l )-XX(I >> 

22s- SLOP E*XXTI1*TYTI1 
PR  INI* (2 X. A.I2.A.F9.4)'* 

1 •THE*»d»  •  TH  PERCENTILE  IS»» 

2<(J/100. t-22  >/ SLOPE 
PRINT* 

GO  TO  410 
CM)  IF 

420  CONT INUE 
410  CONTINUE 


00  430  AK=l,R£p 
K=REP*1-AK 

IFIVYCK) .LT. .991THEN 
60  TO  999 
END  IF 
CONTINUE 

SLOPE=f YY(K*1 )-YY (K>)/(XX(K»I )-XX<K>> 

11-— SLOPt*XX(K).YYiK) 

PRlNTM2KtA.F9.41* , 

I »THE  99TM  PERCENTILE  IS*. 

2I.99-ZZ1/SL0PE 

PRINT* 

PRINT* 

PRINT* 

PRINT* 

CONTINUE 

PRINTM2X«F9.4t4XtF9.4*4XtF9.4l« .CSJ.TSJ.ASJ 
ENO 


•SUBROUTINE  TO  EVALUATE  ENDPOINTS 


SUBROUTINE  ENOPT I XX.YY.REP.NUHI 
INTEGER  REP 

DIMENSION  XXTREPl.VYTREP) 

SLOPE=<YY<21-YY<31 I /( XX ( 2 1 -X X< 3) I 
B= YY 12 1- SLOPE • XX (2 ) 

Vls-B/SLOPE 
IFIVl.LT.O.ITHEN 
VI  =0. 

ENO  IF 
XXU  »  =  VI 

SLOPE: IVY tNUN>-V Y INUM*1 >>/ C XXI  NUN) -XX(NUM»1> 1 

BisVVINUH 1-SLOPE*  XXINUM1 

V2M1.-B11 /SLOPE 

XXtREP 1=V2 

RETURN 
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PROGRAM  3 1  OGGI 

Cl  •»  •  t A  •  ■  A  •  A  #  •  »  •  I  »  III  I  I  M  I  I  l  «  III  I  I  «  i  III  .III  k  •  * 

C*  u.  ...  ..  ..  .4  ...  .  ...  *..*». 

C  THIS  PROGRAM  GENERATES  A  POU'R  STUDY  BETWErN  THE  FOLLOWING. 

C  KS.CRAMlH  VON-MISES.  AIOERSOM* DARLING.  AND  CHI-SQUARE  STA1 

C  5000  HEP  ; 

C  THIS  POWER  STUDY  IS  VALID  FOR  THE  GAMMA  0ISTRI8UTI0N 

C  N  -  SAMPLE  SIZE  =  25 

C  SS1=0  IF  SCALE  PARAMETER  THETA  IS  KNOWN 

C  SS1  =  1  IF  THETA  IS  TO  BE  ESTIMATED 

C  SS2  =  0  IF  SHAPE  ( K I  IS  KNOW 

C  • SS2=l  IF  ( K I  I.S  TO  BE  ESTIMATED 

C  SS-J-0  IF  LOCATION  (C>  IS  KNOWN 

C  - SSS=1  IF  C  IS  TO  BE  ESTIMATED 
C  ■  Cl-INITIAL  ESTIMATE  OF  C  (OR  KNOWN  VALUEI 

C  •  TL=INI T I AL  ESTIMATE  OF  THETA  (DR  KNOWN  VALUEI 

C  A1=INITIAL  ESTIMATE  OF  K  (OR  KNOWN  VALUEI 


C  ...  **■  .*  ....  ...  ...  •»*•»*  »  I  * 

C  •  .  . .  •  a.  ...  ..  .(II  III  II  l|  llll1  I  M  .  I>  ••(  ■  .1.1  .  .  • 


C0MM0N/VALUE/P(100  I 
COMMON/R A  Y/T (1001 

COMMON/M AN A/ A. SSI .SS2 , SS3 , H. C 1 , T 1 , A1 , MR 
DIMENSION  FX ( 60  I tF IX  (60  I  • 

2AAWCVM (50001. AAWAO(50O0l*AAKS(S0O0l 
OOUULE  PRECISION  OSEEOt  T .Cl.Tl.Al 
INTEGER  PP 
OSEEO- 25000. 000 
MR  =0 

RWCVM-O. 

RWKS-O. 

RWAO=0. 

NZERO=0 
NONE sO 

710  CONTINUE 

REAO  •SSl«SS2.SS3fCltTl.Al 
PRINT* •SSllSS2tSS3«Cl  .Tl.Al 
PRINT 
PRINT. 

PRINT 

PRINT* (2Xt A.F5.I »« ,* SHAPE =*,Al 

PRINT*  (2X.AI  •  ,*- — - • 

PRINT. 

PP  si  0 
N=PP 
M  =  N 

00  99  KK=1 *5000 

CALL  GGNLG(0SEE0»N»Q.»1. •>! 

00  719  IKsl.N 
P( IK  I  =  1.  P( IK  1*10. 

719  CONTINUE 

CALL  VSRTAIPtNI 
00  3  IJM.N 
T(IJI=P( IJI 
3  CONTINUE 

CALL  GAMMA(CSJ«TSJ|ASUI 
DO  RB  L=1»N 
Ws(T(LI~CSJI/TSJ 


16  4 
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XX=ASJ 

CALL  HOG AH  <N*XX*  PROB •  ICR  I 

F((L> =PR08 

FIX<LI=FX<LI 

IF  (FIX(L  l.tO.O. ITHEN 

F1XCLI=F  l  X I L •♦•OOOl 

N2CRO=N2CRO«l 

END  IF 

IF  <FIX  (L  I  .'0.  1.1  THEN 
F1XILI=F IXtLI-.ODOl 
NONE  =NONC* 1 
ESO  IF 

GO  CONTINUE 
UCVM=0. 

XN=N 

UAO=0. 

TOP=Q. 

BOT=0. 

00  500  I=1,N 
XI  =1 
RL  =  I 

IF<RL/XN-FIX(I»  .GT.  TOPI  TOP: RL/XN-F I x I l I 
IFIFIXCI l-IRL-ll/XN  .GT.  iOTI  JOT:F I  X (II - CRL-1  l/XN 
HCVM=UCVK*lFXtH-«2.  XI-I.I/12.  XNII  -2 
NA0:UA0«<2.  Xl-I  I  •  C  LOGIFUt  1 1  I  *L3Gt  1-F  IX  (N*l-1 1 1  I 
500  CONTINUE 
OIF=T  OP 

IF (80 T  .GT.  OlFIOIF^BOT 

UXS-DIF 

AAKStKKMOIF 

IFIUKS  .GT.  .27721  RyKS=RUKS»l . 
UCWM=UCtfM*l./(12."XNI 
AAUCVM(KKI=UC VM 

IFtWCVH  .GT.  .139HRVCVH=RVCV4*I. 

MAD=  <“MA  0/ XN I «XN 
AAUAOtKK l=WAO 

IF (UAO  . GT.. 8031  I  RMAO=RUAO»l 
99  CONTINUE 

CALL  VSRIATAAKS, 50001 
CALL  VSR TAIAABCVH,5000 I 
CALL  VSRTAtAAVAOf 50001 
PRINT* 

PRINT'  , •  N2Efl 0=  •  .  N2ER0 
PRINT*  »* NONE:  •  .NONE 
PRINT 

PRINT*  » ' SAMPLE  SI2E  =  *iPP 

PRINT 

PRINT 

PRINT*  *• TOTAL  REJECTION  X  FOR  <-S=  »,RVKS/5000 
PRINT* 

PRINT* «• TOTAL  REJECTION  I  FOR  JCVN:  'tRVCVM/5000 
PRINT 

PRINT* i* TOTAL  REJECTION  X  FOR  JAO:  *,RWA3/5000 

PRINT 

PRINT 

PRINTM2X»F9.4*4K»F9.4.4X*F».4 l» tCSJ.TSJ.ASJ 
ENO 


•SUBROUTINE  GANN*  CALCULATES  RLE  PARAMETERS 


C 

c 
c 
c 

SUBROUTINE  GAMMA1CSU.TSJ.ASJ> 

COMMON  /KAT/TdOO) 

C0MMUN/MAHA/N.SS1»SS2,SS3,M,C1 .Tl.Al.NR 

DOUBLE  PRECISION  T .C » THE T A * ALP HA. OL T • OLA*  AL , OLC ,CE» TH. EN.EN.ELNH 
DOUBLE  PRECISION  EMR ,EI *02 T . OT .D2A.OA »D2C.DC ,ENS,6AM,GMA .GAM I » 6MAI 
DOUBLE  PRECISION  GMAI 2 . DEXP. DABS ,DLOG t SLt DGtOGI ,06 12 »SR . SI .OGAM 
DOUBLE  PRECISION  OGAMI »EL .CSJ. T Sj * ASJ. Cl . Tl • A1 tOSEED 
OIMENSION  CiU00>  .ThETA!  1100  >»  ALPHA!  1100  > 

DIMENSION  OLTIbOI tDLAISOI • AL1 50 > .DLC ! 50 > .CE! 50>tTH!50> 

UI=20 
JH=20 
C!1)=C1 
THETA! 1) =  T1 
ALPHA! l»=AX 
9  EN=N 

EM  =  M 

BG  ELNM=O.DO 

EMR -MR 
MRP=MR  *1 

8 7  NM=N-M*l 

00  88  I- NMiN 
EI=I 

88  ELNH=ELNH*DLOG!£I> 

IF!HR>66.89*109 

109  00  110  I =1.MM 
E I  =  l 

110  £LNM=ELNH-OLOG!£I> 

89  00  63  J= 1.11 00 

IF  !J-l>  66.112.111 

111  UJ=J-1 

IF  1U-JI >6.139.139 
139  IF  ! J/JH-JJ/JH 16.6.117 
117  J2=J-2 

J3=J-3 

IF1SS1 >119.119.1 18 

HR  02T=THETA!JJ 1-2.D0 .THETA! J2 > *T HE T A t J3 > 

OT=THETA! JJ>- THETA! J2) 

IF  !02T> 135.119*135 

1 35  NT=OABS! DT/D2T ) 

GO  10  120 

119  NT  =999999 

120  IF  !SS2>122.122.121 

121  02A= ALPHA! JJ >-2.00* ALPHA! J2 > ♦ALPHA! J3 > 

OA=ALPHA! JJ> -ALPHA! J2> 

IF  i02A113fa.122.13b 

136  NA=OABS( 0  A/0 2 A ) 

GO  TO  123 

122  NA=999999 

123  IF  1SS31125.12S.124 

124  D2C=C( JJ)-2.U0*C!U2>*C! J3> 

0C=C! J J> -C !J2> 

IF  !C!JJdS.0-S-T!l>>14Q«l25«125 


'i. 

i 
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* 

140 

« 

IF  !C(JJ)-5. 0-5)125. 125, 141 

i ", 

141 

If  (020  117. 125. 137 

137 

NC=DABS( 0C/Q2C  1 

CO  TO  126 

■*« 

12% 

NC=499949 

1 

»  126 

NS=2«HIN0(NT,NA*NC) 

IF (NS) 6. 6.142 

142 

IFINS-99 9999 >139.6.6 

138 

ENS-NS 

IF  (SS1> 127.127.128 

127 

THET  A( J) =  THE TA(JJ> 

60  TO  129 

I 

128 

THETAt  J>  =  THET  A!JJ)»(DT ♦  • 2500  •  (■:4S«1. 00  >*D2T)*ENS 

THETA! J)=0MAX1! THET A! J),l.D-4> 

124 

IF  (SS2> 130.130 .131 

130 

ALPHA! J) = ALPHA! JJ) 

CO  TO  132 

131 

ALPHA!  J)  =  ALPHA!  JJ  >  *  (DA»  .2503  *(E9S*1.00)«02A)  *ENS 

ALPHA! J) =DMA  XI (ALPHA! J) .1*0-4) 

132 

IF  ( SS3) 133.133.134 

133 

C( J)=C( JJ) 

i*'. 

134 

CO  TO  112 

C!  J)  =C! JJ)+!OC**25DO*!£NS*1*DO  >*D2C)*£NS 

C! J)=ONAX1!C!J)iO.OQ) 

C! J) =DHI N1 (C(J).T(l)) 

IF  10*00 -EHR)*C! J)-T tl))112.&,6 

|  W 

6 

THETA! J)= THETA (JJ) 

IF  (SSI) 13.13.7 

1 

SlsQ.OQ 

00  8  I=HRP,h 

8 

Sl-Sl*T(l ) -C(JJ> 

IF  (N-H«MR)66, 73.74 

* 

73 

THET  A! J) =S1/!EH. ALPHA! JJ) ) 

CO  TO  13 

1 

74 

6NA=CAM( ALPHA (.M)) 

KS-0 

00  108  K=1.S0 

KKsK-1 

6HAI =CAHI ( !T (H)-C ( JJ> )/THET A ( J) « ALPHA ! JJ) ) 

-'HAI2=GAHI  !!TiHTP)-C(JJ>>/MEM(  J)  .ALPHA!  JJ)  ) 

OLT! K)=-EH» ALPHA (JJ)/THETA(J> *51/ (THETA! J > • *2 >*(EN-EH> • ( * (N)-C(JJ) 

1 ).. ALPHA (JJ) *0£XP! !C (JJ) -T( N ) )/T  HE  TA ( J ) >/ ( THET  A! J) • * ( ALPHA! J J) *  1 *U 

2  0)*!  CHA-GHAI  )  )*EMK« ALPHA! JJ) /THE  TA! J) -EMR  *( T ( NRP )-C ( JJ) ) *  *  AL  3H  A (  J J 

3)*0EXP((C(JJ)-I!HHP))/THETA(J))/!Tr(ETA(J).. ( ALPHA! JJ) *1.00 ) .GHA I 2  ) 

»-f‘ 

101 

TH!K)  =  THE  T  A! J) 

IF  (OLT(K))101 .13.102 

KSsKS-1 

IF  IKS +K >105.103.105 

102 

ks=ks* i 

1F1KS-K) 105*104*105 

■i 

103 

THETA! J) =.5DQ*Trt(K) 

CO  TO  10b 

104 

THET  A(J)=1*50Q.TH!K> 

CO  TO  108 

105 

IF  !OLT!K)*DLT!KK))107, 13.106 

106 

KKsKK-l 

CO  TO  105 

-  ^ 

►“  ; 

/■' 

/ 

* 

* 

$ 

167 
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THETA!  J}  =  TH!K}«OL*TIK>*iTHIKl-TH!KKIl/iOLT!KKI-OI.T!KII 

IF  i DABS! THE TAiJl-THIKl l-l .0-4  >13,13,108 

CONTINUE 

ALPHA! J! = ALPHA! JJ! 

IF  ! SS2I 44,44 , IS 

SL=0.00 

00  16  I-HRP.H 

SL =SL*OLOG(T  !I)-CIJJ)I 

KS=Q 

DO  43  K= 1,50 
KK=K-l 

GMA^GAH! ALPHA ( J) > 

IF  !N-M*MR>66*30«21 

GHAI=GAHU ITiMl-CIUJll/THETM J), ALPHA! J>> 
6MAI2=GAHIi!TIH8P>-C!JJI!/THETA( J), ALPHA! JM 
OG=OGAH! ALPHA! J>  > 

IF !N-H*Hfl  >66,77,32 

DLA!K1=-EH*UL06! THETA! J> ) * SL-E 4* 06/6MA 
GO  TO  7B 

OGI-OGAM I!!! !N> -C 1 JJ11 /THETA ( J), ALPHA! J> > 

DGI2-DGAHI  i!TIMRP)-CI JJ )>/T8ETA!J),ALPHA!J>> 

DLA(K>--EH*DLOG!THETA!J)>*SL-EN* OG/GHA*! EN-EH} • 1 OG-DG l 1 / i GHA-GH AI 1 
1 *EHR»DL0G ! THETA!  J) l*EHR *DGI2 /64 AI 2 
AL  !K  1  =  ALPHA!  <11 
IF  ! OLA! K >134 ,44 ,40 
KS=KS-1 

IF  !KS«K 170,41 ,70 
KS=KS*1 

IF  IKS-KI70, 42,70 
ALPHA! J>=. 500* AL1K1 
GO  TO  43 

ALPHA! U>=l«5O0*AL!Kl 
GO  TO  43 

IF  !OLAIK1*OLA!KK1 172,44,71 

KK=KK-l 

GO  TO  70 

ALPHA! J13ALIKI*DLA!KI«!AL!KI-AL!KK11/IDLAIKK1-0LA!K1I 

IF  IOABSiALPHAIJI-ALiKll-l.O-4144,44,43 

CONTINUE 

C!UI=C! JJ1 

IF  ISS31112, 112*45 

IF  ! 1 .00 -ALPHA! J >  >79, 143, 143 

IF  !SS1*SS2>57,57,79 

IF  IN>N)66,b 3,46 

GNA =G AN!  ALPHA! J1 1 

KS=0 

00  56  K=1,S0 

KK=K-1 

SR SO. 00 

00  69  I:NRP,M 

SR=SR*1.00/!T! Il-Cl  Jll 

IF  !N-H*NR>66,80,ai 

OLC! K>=! 1.00 -ALPHA !U>)*SR*EH/ THE T A! J> 

GO  TO  62 

GNA I =G AH  I! !T !H>-C! J>> /THETA!  J>,  ALPHA! Jll 
GNAI2  =  GAHI1! T ! HRP > -C I J > > / T HE T A  I J > (ALPHA! Jll 

OLCIK1M 1.D0-ALPHA! Jl 1 *SR*IEN-iR8 1/THETA! J) *! EN-ENI* IT INI-CIJI1**! 
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« 

1  ALPHA! J) -1 .00 )*DEXP!-!T !N)-CiJ))/THCTAiJ>)/( THET  A !J) • ‘ALPHA! J)  •  <  G** 
2A-GMAI) >-EMH*IT<MRP>-CCJ>>**! ALPHA! J>-1.00>*OEXPC-<T IMRPI-CI J) )/TH 
3ETAIJ))/! THE  TAIJ) •* ALPHA! J) *64 At  2) 

82  C£IK)=CIJ> 

51  IF  IDLCIKI)9Q*112,91 

90  XS=KC-1 

IF  I KS*K ) 54*52  *54 

91  KS=KS*1 

IF  IKS— K 154*53*54 

52  Cl«l)  =«5D0  >CE  (K  > 

GO  TO  6a 

53  Ct  J>=CEiK>«.5D0MTU>'CEIK>> 

GO  TO  68 

54  IF  IOLC!K)‘OLC!KK 1)67 *112*55 

55  KK-KH-l 
GO  TO  54 

67  Cl J)=CE!K>*OLCtK>  *ICE!K)-CEf  KO>/IOLCtKK>-DLCtK>> 

68  IF  I0ABS!C(J>-CCIK>>-1. 0-4)112, 112*56 

56  CONTINUE 
GO  TO  112 

57  CIJ)=TI1) 

112  IF  INR>66*113*S8 

113  DO  115  1=1, M 

IF  ICIJ)«1.0-4-T!I 1)116*114*114 

114  HR=MR*1 

115  CI1)=TI1) 

116  IF  INR >66,58*86 

58  Sl=Q.OQ 
SL=0.D0 

00  92  I  =  MRP  *  M 
S1=S1*T( I  l-CIJ) 

92  SL=SL«OLOG(III)-CIJ)> 

GMA=GAMI ALPHA! J) ) 

IFIN-M*MR >66*98*96 

96  GMAI=GAMI 1 1 T  IN  >-C  !.)>  > /THET A!  ,1)  ,  ALPHA  I J)> 

GNAI2=GAMIII TIHRP)-C!J)>/THETAIJ>,ALPHAIJ>> 

98  CL=ELNM-EM*OLOGI GM A) -EM •ALPHA! J)‘ DLOGITHETAtJ))*! ALPHA IJl-l* DO >‘SL 
1 -Sl/THEI A I J) 

IF  IN-H«M*)66, 100*99 

99  EL=£L+ IEN-EN ) • tOLOCIGMA-GNAI ) -DLOGIGNA) ) 

1  ♦EM8  •  ALPHA!  JXOLUGI  THETA!  J) > *ENft *0LOGI GNAI2) 

100  CSJ=CIJ) 

TS«1=  THET  A I  J) 

A$U= ALPHA (J) 

IF  I J-2) 63 *60 ,60 

60  IF  I0ABSIC!J)-CIUU))-1. 0-4)61*61*63 

61  IF  lOABSITHETAIJI-THETAIJjn-l. 0-4)62*62*63 

62  IF  IOABSI ALPHA! J) -  ALPHA! JU) )-l»D *4 14*4*63 

63  CONTINUE 

4  CONTINUE 

66  RETURN 

ENO 


If'.; 


FUNCTION  SAN 


74/74  OP 7=0 


FTN  5.1*524 


C  •.••*••••..••*••»*••*••*•••••**•**•*•*••*•••••*»•••••■••••■•■•■■ 

C  *•• . . . . . . . . . 

DOUttlC  PRECISION  FUNCTION  GAM  f  I 
DOUBLE  PRECISION  6*2*OLOG»DEXP,t 
Z-1 

6=0.00 

1  IF  (2-9. 00)2,2,3 

2  G=G-DLOG(Z> 

2=2*1. DO 

60  TO  1 

3  GAH=G*(Z-.S001*DLOG(Z!-Z*.500*)t.  0G(2. 00  *  3.14  1  592  6  5  3589T9  900  7*1.00/ 
l (12. 00*2 1 -1.00/( 360.00 *2* * 31 *1.00/ (1260.00*2* *51-1.00 /(l 680.00 *2* • 
2  71*1. 00/ ( 11 bo .00*2**91-61 1. 00/(360360.00 *2* • 1 11*1.00/ (156.00*2**13 
31 

6AH=0£XP(GAM1 

RETURN 

EDO 
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DOUBLE  PRECISION  FUNCTION  0644(7! 

DOUBLE  PRECISION  06*2*T«DL0S*SAN 
2=7 

06=0.00 

1  IF  (2-9.0012*2*3 

2  06=06-1.00/2 

2=2*1.00 

60  TO  1 

3  06AH=OG*(2-.500!/2*DLOG(2!-1 .00- 1 .00/(1 2.00 *2**21* 1 .00/( 1 20 .00 *7* • 

1  41-1. 00/(252. DO *2**6!*1. 31/(240. 00* 2* .0)-l. 00/(1 32. DO *2**  10 1 

2  *69 1.00/(32 760.00*2* *12 1-1. 00/(12. DO *2* *141 
064N=06AN*6AN(T1 

RETURN 

CNO 


FUNCTION  0G4NI  74/74  0PT=0 


FTN  5.1*524 


1 


OOUBlE  PRECISION  FUNCTION  0SANl(N»21 

DOUBLE  PRECISION  U»7»M*2»SU»Ell 

DI MCNSI ON  U( 501*7(501 

U(ll=M**2*DLOG(W»/2 

7(  1 1  =N  **  2/2  •  *2 

SU=U(1 1-7(11 

DO  1  L=2»50 

LL=L-1 


U(L1=(”U(LL1 • M/ELL) *(2*ELL 

7(L!=-7(LL»*7*(/*CLL-1  .001 

SU=SU*U(L1-7(L1 

OGANl-'SU 

RETURN 

CNO 


-1.091/(2«CLLl 

•*2/((2*ELLI**2*ELLl 


17(i 


FUNCTION  6AHI 


14/74  0PT=0 


FTN  S«l*524 


DOUBLE  PRECISION  FUNCTION  SAMI CU»2) 
DOUBLE  PRECISION  U«U*2«SU»ELL 
DIMENSION  U<50> 

U(l»=U**2/2 

SU=Utl» 

DO  1  L=2 1 50 

LL=L-1 

£LL=LL 

U«L)=<-UCLL>/ELLI«U»(2*ELL-l.Qll/I2*ELL> 
I  SU=SU*l/I  L  ) 

GAHI =SU 
RETURN 

END 
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